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Duality relations for the ASEP conditioned on a 
low current 
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Abstract We consider the asymmetric simple exclusion process (ASEP) on a finite 
lattice with periodic boundary conditions, conditioned to carry an atypically low cur¬ 
rent. For an infinite discrete set of currents, parametrized by the driving strength sk, 
K > l,we prove duality relations which arise from the quantum algebra t/^[g[(2)] 
symmetry of the generator of the process with reflecting boundary conditions. Using 
these duality relations we prove on microscopic level a travelling-wave property of 
the conditioned process for a family of shock-antishock measures for N > K parti¬ 
cles: If the initial measure is a member of this family with K microscopic shocks at 
positions (.r:i,... ,xk), then the measure at any time f > 0 of the process with driv¬ 
ing strength sk is a convex combination of such measures with shocks at positions 
(yi,... ,yK)- which can be expressed in terms of A'-particle transition probabilities 
of the conditioned ASEP with driving strength Sfj- 


1 Introduction 

In the asymmetric simple exclusion process (ASEP) Il3^l2^l27ll?4l each lattice site 
k on a lattice A = (1,... ,L) is occupied by at most one particle, indicated by occu¬ 
pation numbers ri{k) €S = {0,1}. We denote by T] = ( 77 ( 1 ),..., T 7 (L)) € Q = a 
configuration 77 of the particle system. Informally speaking, in one dimension parti¬ 
cles try to jump to the right with rate r = wq and to the left with rate £ = . The 

jump attempt is successful if the target site is empty, otherwise the jump attempt 
is rejected. The invariant measures of the ASEP with periodic boundary conditions 
are well-known: For fixed particle number N these are the uniform measures. From 
these one can construct the grandcanonical Bernoulli product measures with fugac- 
ity z = p/(l — p) where p = N/L is the particle density on the torus. For these 
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measures, where each lattice site k is occupied with probability p independently of 
all other sites, one has a stationary particle current j* = {r — i)p{l— p), correspond¬ 
ing to an expected mean time-integrated current (7(f) )/f = j*- 

In the context of macroscopic fluctuation theory 0 one is interested in condition¬ 
ing the process on fluctuations around some atypical mean time-integrated current 
j 7 ^ j* ■ A question of fundamental interest is then which macroscopic density pro¬ 
file is most likely to realize such a large deviation of the current inside a very large 
(more precisely: infinite) time interval of conditioning. This large-deviation problem 
thus concerns an untypical ensemble of trajectories of the process. This ensemble is 
usually not defined by j but via Legendre transformation in terms of the canonically 
conjugate driving strength s{j) with s{j*) = 0. Interestingly, for conditioning on a 
lower-than-typical current (i.e., for s < 0), it was found by Bodineau and Derrida ||9l 
for the weakly asymmetric simple exclusion process that there is a dynamical phase 
transition: For currents slightly below the typical value j* the optimal macroscopic 
profile is constant as it is for j*. However, below a critical threshold jc < j* (corre¬ 
sponding to some Sc < 0) the optimal macroscopic profile is a travelling wave with 
a shape resembling a smoothened shock/antishock pair. 

More recently, in a similar setting, but for finite duration t of conditioning, the 
microscopic structure of a travelling wave in the ASEP (not weakly!) was elucidated 
in detail for a specific choice of negative driving strength 0 : One considers a certain 
family of inhomogeneous product measures indexed by a lattice site k, where the 
microscopic density profile as function of the position on the lattice has a density- 
jump at position k on the torus, analogous to a shock on macroscopic scale. At time 
t = 0 N particles (N arbitrary) are distributed according to the restricted measure 
l^ic^f.ri{k),N- Then at any future time f > 0 of the conditioned dynamics the 
measure is a convex combination pj^ (f) = of such measures. The 

weights c{l,t\k,Q) are the transition probabilities of a single biased random walk, 
thus suggesting that a shock in a macroscopic travelling wave performs a biased 
random walk on microscopic level. 

In this work we trace back the mathematical origin of this rigorous result to 
certain algebraic properties of the generator of the process. Then, using these prop¬ 
erties, we go beyond a to derive a family of duality relations that allow us to con¬ 
struct more complex microscopic structures corresponding to more general macro¬ 
scopic optimal profiles. The starting point is the well-known fact that for reflecting 
boundaries, where the process is reversible, the generator of the process commutes 
with the generators of the quantum algebra t/^[g[( 2 )] ([T]. This fact has been used in 
ED to construct the canonical reversible measures and in to derive self-duality 
relations for the unconditioned ASEP. 0 

On the torus, however, the symmetry of the generator under t/^[g[(2)] breaks 
down. Nevertheless, some time ago Pasquier and Saleur ll^ found intertwining 
relations involving the generators of Uq[Q\.{2)\ and the Heisenberg quantum Hamil- 

* We mention that the deep link between duality of Markov processes and symmetries of its gen¬ 
erator, first noted in (32), that we exploit here was given a systematic abstract treatment in ISD. 
More recently many concrete symmetry-based dualities for interacting paificle systems were de¬ 
rived using this approach (29][II][l3][Il][l0][l2l|7]|25][l5l. 
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tonian with a boundary twist. This quantum Hamiltonionan operator became later 
to be known to be closely related to the generator of the conditioned ASEP ifTSll . 
Here we present a new proof for the results of ll^ (and correct some typos there) 
for reflecting boundaries and make use of intertwining relations of Go) (correcting 
some typos also in that paper) to derive an infinite discrete family of duality rela¬ 
tions for the ASEP with periodic boundary conditions. These new duality relations 
apply to the process conditioned on fluctuations around some untypically low mean 
time-integrated current. 

The simplest of these duality relations proves that the homogeneous Bernoulli 
measure is the invariant measure for the unconditioned ASEP with periodic bound¬ 
ary conditions. The derivation of this well-known fact from the t/^[g[(2)]-symmetry 
of the process with reflecting boundary questions is remarkable in so far as it raises 
the interesting question whether one can construct the matrix product measures 
II 20 I [191 of the periodic multi-species ASEP from the L(j[gt(«)]-symmetry of that 
process with reflecting boundaries|[T]|6l. 

Prom the non-trivial higher order duality relations we obtain an infinite discrete 
family of new microscopic “travelling waves” for the conditioned process. 


2 Definitions and notation 

It is convenient to work with the quantum Hamiltonian formalism ll28l |34l where 
the generator of the process is represented by a matrix which in a judiciously cho¬ 
sen basis turns out to be closely related to the Hamiltonian operator of a physical 
quantum system. We first introduce some notation and then describe in some detail 
the tools required for the quantum Hamiltonian formalism for the benefit of readers 
not familiar with this approach. 


2.1 State space and configurations 

We say that a site k € A is occupied by a particle if t] (k) = 1 or that it is empty if 
rj (k) = 0. The fact that a site can be occupied by at most one particle is the exclusion 
principle. Occasionally we denote configurations with a fixed number of N particles 
by rjN- The set of all configurations with N particles is denoted Qn- We also define 

t;(k):=l-77(k) (1) 


and the particle numbers 

L L 

Niri)=Y,ri{k), V{ri)=Y,v{k)=L-N. 

k=l k=\ 


(2) 
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A useful alternative way of presenting uniquely of configuration rjN is obtained 
by labelling the particles consecutively from left to right (clockwise) by 1 to A and 
their positions on A by x, mod L. A configuration 77 is then represented by the set 
X : = {x : 77(x) = 1}. We call this notation the position representation. We shall use 
interchangeably the arguments 77 , x, for functions of the configurations. When the 
argument is clear from context it may be omitted. We note the trivial, but frequently 
used identities A( 77 ) =A(x) = |x| and 


7V(x) 

77(k)= ^ 5,,,,. (3) 

!=1 

For a configuration 77 = x we also define the number Ajt (77 ) of A-particles to the left 
of a particle at site k 

k-\ N{ri)k-\ 

Nk{ri) ■■= E ^7(0 = L L (4) 

!=1 i=l /=! 

Furthermore, for 1 < k < L — 1 we define the local permutation 
7r'^<^+i(77) = {77(l),...77(k-l),77(k + l),77(k),77(k + 2),...,77(L)}=: 77 ^^+', (5) 
and for k = Lwe define 

77^1(7?) = {77(L),...,77(k),..., 77 ( 1 )}=: 77^^ ( 6 ) 

The space reflection is defined by 

/?(77) = {77(L),77(L-1),...,77(1)} (7) 

corresponding to R{ri{k)) = rj{L +1 — k) for the occupation numbers. 


2.2 Definition of the ASEP 

For functions / : —7 C the ASEP 77 , with periodic boundary conditions and hop¬ 

ping asymmetry q is defined by the generator 

■^/(j?) in')-fin)] (8) 

where the transition rates between configurations 

w (77 -4 ri') = Y.w’^^+\ri)8^, ^kk+i (9) 

k=l 


are defined in terms of the local hopping rates 
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=w[qr]{k)v{k+l)+q 't)(^)77(^+1)] . (10) 

The prime at the summation symbol ® indicates the absence of the term Vj' = Vj 
which is omitted since w(77 rj) is not defined^ The transition rates are non-zero 
only for a transition from a configuration 77 to a configuration 77 ' = defined 

by (01. 

We shall assume partially asymmetric hopping q ^ 0,1,°°. The constant w ^ 
0 sets the time scale of the process. On the torus we identify increasing order of 
the lattice index with the clockwise direction. In the case of reflecting boundary 
conditions no jumps from site 1 to the left and no jumps from site L to the right are 
allowed. Increasing order of the lattice index is identified with the direction left to 
right. The upper summation limit L in (0 has to be replaced by L — 1, giving rise to 
a generator that we denote by .if. 

In order study fluctuations around some untypical integrated current, parameter¬ 
ized in terms of the driving strength s, we define the weighted transition rates 

= w[^e*77(k)7;(k-f I)-f^-^e^"7;(k)77(k-f I)], I <k<L-l(lI) 

This leads us to define the weighted generators 

:= Ewf+i(77)/(77“+')-w“+H77)/(77) (13) 

yt=l 

The weighted generators give a weight e^ (e^®) to each particle jump to the right 
(left) anywhere on the lattice and for the process with periodic boundary conditions 
an extra weight e® (e^®) to each particle jump to the right (left) across bond (L, 1). 
Thus each random trajectory of the process is given a weight where 7(f) 

is the time-integrated total current, i.e., the total number of all particle jumps to 
the right up to time f minus the total number of all particle jumps to the left up to 
time f and 7^(7) is the time-integrated current across bond (L, 1). Notice that the 
diagonal part of the weighted generator does not depend on the driving strength 
s or s, reflecting the fact that the random times after which jumps occur remain 
unchanged. For details on this construction see e.g. II 221 123 [161 specifically for 
the present context |[35l. 

We fix more notation and summarize some well-known basic facts from the the¬ 
ory of Markov processes. For a probability distribution P{rj) we denote the expecta¬ 
tion of a continuous function f{rj) by {f)p := Y.ri Th® transposed gen¬ 
erator is deflnedby f{T]) L^^/g§z./(t 70 -^l 7 )'(fl) where With 

this definition ([8]) yields for a probability distribution P{t]) the master equation 


^ When the summation is over we shall usually omit the set under the summation 

symbol and simply write . 
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■^^P{ri) = ^ n)P{ri')-Hri ^ i?')^(i7)]- (15) 

An invariant measure is denoted 7t*{ri) and defined by 

^'^n*{ri) = Q (16) 

and the normalization Ltj ^*(^) = 1 - unnormalized measure with the property 
(IThl) is denoted n{rj). The time-reversed process is defined by 

:= EV-(J? ^ J7')[/(^7') -/(^7)] (17) 

7)' 

with —>■ T]') = w( 77 ' ri)7t{ri')/7z{ri). The process is reversible if = 

jSf which means that the rates satisfy the detailed balance condition n{rj)w{ri —> 
rj') = w(ri' ri)7z{ri'). A probability distribution satisfying the detailed balance 
condition is a reversible measure. It is easily verified that the ASEP with reflecting 
boundary conditions is reversible with reversible measure 

;r( 77 ) = 


for any ^ S R. 

We define the transition matrix H of the process by the matrix elements 


_ r-w(77 -> 77') 77 ^77' 

{'L'nMPri') ri = ri'. 

with w{ri —!■ 77 ') given by (|9]l. One has 

r»' r»' 


(19) 


( 20 ) 


Notice that here the sum includes the term 77 ' = 77 . In slight abuse of language 
we shall also call H the generator of the process. Analogously we also define the 
weighted transition matrix where the off-diagonal elements are replaced by the 
weighted rates (dUl, (Ell. 

For an unnormalized stationary distribution we define the diagonal matrix n with 
the stationary weights 7r(ri) on the diagonal. For ergodic processes with finite state 
space one has 0 < n(rj) < for all 77 . In terms of this diagonal matrix we can 
write the generator of the reversed dynamics as The reversibility 

condition = H then reads 


= ( 21 ) 

Therefore, if one finds a diagonal matrix with the property (l2TI) then this matrix 
defines a reversible measure. 
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In order to write the matrix H explicitly we assign to each configuration rj a canon¬ 
ical basis vector 1 77 ). We choose the binary ordering 1 ( 77 ) = 1 +Y!k=i T]{k)2^^^ of 
the basis. Defining single-site basis vectors of dimension 2 

| 0 ):=(J). | 1 ):=(") ( 22 ) 

one then has 1 77 ) = |77(1)) (g)... (g) |77(L)) where (g denotes the tensor product. 
These basis vectors span the complex vector space of dimension d = 2^. We 

also define transposed basis vectors ( 77 1 := 177 and the inner product (v' | w) := 
Ir, v{ri)w{ri). 

Furthermore we define the two-by-two Pauli matrices 


CJ' := 





(23) 


and the two-dimensional unit matrix 1. From these we construct 

C7± == n = b = i(l-bcT'). (24) 

These matrices satisfy the following relations: 

= < 7 ^( 7 + nV=Q, Vn = Q 

(j+u = a+, 7!(7+ = 0, a+f) = 0, '0a+=:C7+, (25) 

a^n = Q, ha^ = a^, a^v = va^ =0. 

With the occupation variables ([T]i for a single site we have the projector property 


7i|77) = 77|77), 7)|77) = 77|77). 


(26) 


Having in mind the action of these operators to the right on a column vector, we 
call O' a creation operator, and annihilation operators. When acting to the left 
on a bra-vector the roles are interchanged: (J+ acts as creation operator and o as 
annihilation operator. 

For L > 1 and any linear combination u of these matrices we define the tensor 
operators mj. := 1 ®^^' (g m (g By convention the zero’th tensor power of any 

matrix is the c-number 1 and m®' = u. We note that also the tensor occupation 
operators act as projectors 


N(ti) 

>hl V ) = v(k)\v ) = E ^-^.-41 v), (27) 

!=1 

with the occupation variables 77 (^) or particle coordinates v, respectively under¬ 
stood as functions of rj. The proof is trivial: The first equality is inherited from (|26] | 
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by multilinearity of the tensor product, the second equality follows from ([3]). Mul¬ 
tilinearity of the tensor product also yields uj^Vk+i = 0 [(« 0 1)(1 0 v)] 0 

and the commutator property Uj^vi = viiik 
for k^l. For k = I one has relations analogous to (l25l l. 

It turns out to be convenient to introduce parameters a = and j3 = e' and 
express for periodic boundary conditions the weighted generator as H{q,a,P) with 
the convention H{q,q, 1) = H for the unweighted generator. Similarly one writes 
H{q,a) for the weighted generator with reflecting boundary conditions with the 
convention H{q,q) = H. With these definitions the weighted generators H{q,a) 
and H{q, a,j5) defined by (fOl l and (fl4l i resp. become 

L-l 

Hiq,a) = Y,hk,k+i{q,a) (28) 

k=l 

H{q,a,l5) = H{q,a) + hL,i{q,ci,l5) (29) 

with the hopping matrices 




hL,l{q,a,|i) = -w[apa+a|^^^-qn|,Vk+l + iap) ‘OA+i)] • 

(31) 

It is useful to introduce the space-reflection operator R defined by 


RukR ^ (32) 

for local one-site operators ut and the diagonal transformations 

y( 7 ) = ^ ^33^ 

W{z) = (34) 

with the number operator iV = Y!k=i ^k- We note the properties 

RV{y)R-^ =V-\r)=Vir^), (35) 

RW{z)R-' = Wiz)=W-'iz-^), (36) 

RH{q,a,l5)R-^ =Hiq,a-\li-^) (37) 

H^{q,a,l5) = H{q,a-\l5-^) (38) 

WH{q,a,l5)W-' =H{q,a,l5). (39) 

Moreover, the transformation property 

V{Y)(yj^V-\Y) = 7±i(2<:-L-!)(y± (40) 


yields 
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V{Y)H{q,a,l5)V '( 7 ) =//(^,a 7 ^j3/) (41) 

V{y)H{q,a)V-\j) = H {q,ay-^). (42) 

Thus for periodic boundary conditions global conditioning and local conditioning 
are related by a similarity transformation, while for reflecting boundary conditions 
the conditioning can be completely absorbed into a similarity transformation. One 
also finds with y = q^, a = q the reversibility relation V{q^)H{q)V^\q^) = . 

By (ISTT i this shows that ft = V^^{q^) is the matrix form of the reversible measure 

® with^ = -L-i nail 

For driving strength sq := — In (q) corresponding to a = 1 one finds for 1 <k < 
L-1 


hk.k+i 1) = -^ (T^+i + C7^’0-^+1 + 4 -1) + h{G^^ - )] (43) 

with 

^ = h = (44) 

and the unit-matrix 1 of dimension 2^. Notice that for periodic boundary conditions 
the local divergence terms — (7|^j cancel. For reflecting boundaries the local di¬ 
vergence term contributes opposite boundary fields h{ai — dj"). With the further 
choice So = 0 corresponding to )3 = 1 the weighted generator ( |29] | becomes the 
Hamiltonian operator1,1) of the ferromagnetic Heisenberg spin-1/2 quantum 
chain, while for )3 7 ^ 1 one has the Heisenberg chain H{q^ l,j3) with twisted bound¬ 
ary conditions . 

Since particle number is conserved the process is trivially reducible. For each 
particle number N one has an irreducible process t on the state space We 
define the projector 

Ia'— L l’7)(J7l (45) 

rieOi^/ 

where we have used the quantum mechanical ket-bra convention | rj) (rj | = | ® 

(77 I for the tensor product of two vectors. Thus one obtains the generator 

HN{q,a,P) :=iN H{q,a,l3) In (46) 

for the N-particle weighted ASEP. 

Notice that 1^ acts as unit matrix on the irreducible subspace corresponding to 
particle number N. The unit matrix 1 in the full space has the useful representation 

1=L|J?)(J?I- (47) 

Tien 


^ This is equivalent to Eq. (2.14) in (33), which, however, has a sign error and should read = 

y-W. 
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2.4 The quantum algebra f/^[0((2)] 


The quantum algebra Uq[Q[{2)\ is the q'-deformed universal enveloping algebra of 
the Lie algebra g[(2). This associative algebra over C is generated by i = 1,2 
and with the relations MM 


[L,-,L,]=0 

q-q ‘ 

Notice the replacement q^ ^ q that we made in the definitions of CD. 

It is convenient to work also with the subalgebra Uq[s\{ 2 )]. We introduce the 
generators N and V via = Li, = L 2 and define 

S^ = i(N-V) (51) 


(48) 

(49) 

(50) 


and the identity I. Then the quantum algebra Uq[5i{2)] is the subalgebra generated 
by and S* with relations 




[s+,s-] = : 




■q 


-2S= 


q-q- 


(52) 

(53) 

(54) 


Observing that N + V belongs to the center of t/^[0[(2)] ll24l one sees that t/^[sl(2)] 
is a subalgebra of t/^[gl(2)]. 

It is trivial to verify that t/g[gl(2)] has the two-dimensional fundamental rep¬ 
resentation S* —> (7^, N —?> «, V 0 given by the matrices (l24li . Then (7^ and 
(7^/2 form the two-dimensional fundamental representation of t/^[sl(2)]. Reducible 
higher-dimensional representations can be constructed using the coproduct ll24l 


4 (S±) = (g) q-^" + q^" ® (55) 

Zi(S^) = S"(8)l-fl(8)S^ (56) 


By repeatedly applying the coproduct to the fundamental representation we obtain 

S^{k) = ^2 ‘27-2 Tk+i (J± 


One has 


(57) 

(58) 
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r q^^S^{l)S^{k) 

S^{k)S^{l) = I 0 

[q^^S^il)S^ik) 


(59) 


Thus the spatial order in which particles are created (or annihilated) by applying 
the operators {k) gives rise to combinatorial issues when building many-particle 
configurations from the reference state corresponding to the empty lattice. 

From the coproduct one obtains the tensor representations of C/^[s[(2)], denoted 
by capital letters, 

L L 

S^='^S\k). (60) 

k=\ k=\ 

For the full quantum algebra L(j[gl(2)] the tensor generators are and N = rik, 
V = Vk- The unit I is represented by the 2^-dimensional unit matrix 1 := 1®^. 

For reflecting boundary conditions the Heisenberg Hamiltonian H{q, 1) is sym¬ 
metric under the action of L(j[g[(2)] (Ull^. This symmetry property is the origin of 
the duality relations derived in ll^ and will also be used extensively below. In fact, 
for 1 < k < L— 1 one has [hk,k+i (?; 1) j 5^] = [hk,k+\{<1, Ct) j S^] = Oj which imply 

[H(^,l),5±] = 0 (61) 


and, equivalently to (l39l) . the diagonal symmetries [H{q, a), iV] = [H{q, a), V ] = 
0 , thus giving rise to the t/^[g[( 2 )] symmetry of H{q, 1 ). 

We stress that [liL i(q^ 1), 5^] 7 ^ 0. One the other hand [/li , 5‘] = 0. 

Hence for periodic boundary conditions the symmetry breaks down to only a resid¬ 
ual t/(l) symmetry [H{q,a,l5), S'] = Q generated by S', which corresponds to par¬ 
ticle number conservation since the z-component of the total spin S' is related to the 
particle number operator N through S' = L/2 — N. 

We also define 

L 

S^{q,a) = Y,Sf{q,a) (62) 

k=l 


where 


sUq,ci) = 




(63) 


The diagonal transformation (l40l) and the defining relation (l53l l yield 


V{Y)S^{q,a)V \Y) = S^{q,aY (64) 

W(z)S^{q,a)W-'{z) =z'^'S^{q,a). (65) 


Notice that 1) = S^ as defined in ( l60t . Hence S^{q,a) and S^ also form a 
representation of Uq\5\.{2)\. Since according to (IhTT i H{q, 1) commutes with the gen¬ 
erators 5^ = S^{q, 1) we conclude from (l42li that//(o', a) commutes with S^{q, a), 
which together with S^ form an equivalent representation of L(j[sl(2)]. In particular, 
the generator of the ASEP with reflecting boundary conditions H = H{q,q) com¬ 
mutes with S^ := S'^{q,q). 

We note that 
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{s^iq,a))^ = S^{q,a ') (66) 

(67) 

To prove the second equality one uses RS^{q,a)R^^ = ,a^^) which 

comes from (l3^ . 

Finally we introduce the symmetric ^-number 


r-q 

q-q- 


( 68 ) 


for q,q^^ ^0 and x S C. This definition can be applied straightforwardly to finite¬ 
dimensional matrices through the Taylor expansion of the exponential. For integers 
we also define the ^-factorial 


],! 


1 n = 0 

UU[k], n>l. 


(69) 


2.5 Duality in the quantum Hamiltonian formalism 

For self-containedness we briefly review how to express expectation values using 
the matrix representation of the generator which allows to state the notion of duality 
in a neat matrix form llJTl l2Tll . 

A probability measure P(r]) is represented by the column vector 

= (70) 

v 

Next we define the summation vector 

(^h=E(j?l (71) 

n 

which is the row vector where all components are equal to 1. The expectation {f)p 
of a function /(t]) with respect to a probability distribution P{ri)is the inner product 

{f)p = {f\P) = {^\f\P) (72) 

where 

/:=E/(’?)I’?)(’7I (73) 

v 

is a diagonal matrix with diagonal elements f(rj). Notice that 

f(v) = (vlflv} = (slflv)- 


(74) 
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One obtains the diagonal matrix / ( |7^ corresponding to a function f{rj) by substi¬ 
tuting in f{ri) the variable rj (k) by the diagonal matrix 

For a Markov process rjt the master equation (flST l for a probability measure 
P{r]t) :=Prob[Tj, = tj)] reads 


A|P(f)) = -//|P(f)) (75) 

which implies 

\P{t))=e-^‘\Po) (76) 

for an initial probability measure 7b(i?) = Pi'Ho) at time f = 0. We write the expec¬ 
tation of a function /(?];) as 

(/(O) :=L/(J7)^(j?r) = L/(J?)(J7 le-'^'lFo) = |/V"'|Fo)- (77) 

n v 

If the initial distribution needs to be specified we use an upper index (/(f) )^. Nor¬ 
malization implies (i |P(f)) = 1 for all f > 0 and therefore (s |// = 0. A stationary 
distribution, denoted by | tt* ), is a right eigenvector of P[ with eigenvalue 0, i.e., 
H\n*) =0 and normalization (s|fr*) = 1. For the ergodic subspaces with fixed 
particle number A it is unique. 

In order to introduce duality we consider a process i*, with generator PI and a 
process x, with generator G which may have different countable state spaces £2^ and 
£2b- Consider also a family of functions > C indexed by x S and a family 

of functions g'’ : Qg C indexed by ^ such that/'/^) =g^ (x) =: D(x,^). Let 
the process start at some fixed ^ € £2^ and let X; start at some fixed x S £2b. Then 
the two processes are said to be dual with respect to the duality function D{x,^) if 

ca 

= (78) 

As pointed out in m] this property can be stated neatly in terms of the generators 
as 

DH = G^D (79) 

where the duality matrix D is defined by 

77= E E D(x,^)|x)(^| (80) 


By construction one has D{x,^) = (x |Z)| ^ ). 


2.6 Shock/Antishock measures 

In vector notation a product measure with marginals p*. is a tensor vector | {p,}) = 
|pi) 0 ... ® |pi) with the single-site column vectors jp*.) = (1 — Pk^PkY■ It is con- 
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venient to introduce the local fugacity 


Pk 




i~Pk 


( 81 ) 


and write the product measure in the form | {z;}) = |zi) ® ... 0 \zl)I'^l with |z^) = 
(1,ZA:)^ and normalizationZ/, = HL \[^k/{}+Zk)\- 

Specifically, we define for a set x of lattice sites with cardinality K = |x| the 
following family | Vx) = |zi) Zi ■ ■ ■ (8i \zl)/2l of shock/antishock measures, or SAM 
for short, in terms of the fugacities 

j zq^’^ for xi<k<xt+i, / G {0,1,1} 
for/tGx 


with XQ = 0 and xk+i = L + 1. On coarse-grained scale with ^ = k/L and = 
XijL the macroscopic density profile p(^) corresponding to the fugacities Zk has 
discontinuities at with constant fugacity ratios z^ jzj = where zf = z^r'±i ■ 

In forward (clockwise) direction and for ^ > 1 this is an upward step, corresponding 
to a shock profile for the ASEP (with positive bias ^ > 1). Between site L and 
site 1 there is downward jump with fugacity ratio q^^^. On macroscopic scale this 
constitutes an antishock at position = — (1 + K')/2 mod 1, hence the term 

SAM. These shock measures are closely related to the shock measures defined in 
||4l and also to the infinite-volume shock measures studied in fSl where the shock 
positions x, are occupied by second-class particles. 

With a different normalization factor the general SAM (l82l i with constant fugac¬ 
ity jumps q^ can be written as 




K 

) =n^ 

f=i 




|Vx). 


(83) 


Here 

|z):=k)®^ (84) 

for fG = 0 is the unnonnalized homogeneous product measure corresponding to the 
Bernoulli product measure |p) = |z)/(l-l-z)^ where p = z/ ( 1 -f z). 

From the SAM defined by dS^ we construct a second type of SAM’S using the 
transformations (l33l l and (l34l i 




\p.) ■.= Z-’^V{q=L)Y\W{q^)\P,) 


K 

n 

y=i 


-\l i:f=i"<+ eL +1 «,■ 

Z ({ J i 


l^)- 


We illustrate the definition for K = I and K = 2. 
For K = I the SAM (l85l l reduces to 


(85) 

( 86 ) 
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with 1 <x<L. This corresponds to local fugacities 


Zk 


l(k-x)+L 

zq z for 1 <k <x 

oo for k = x 

2{k-x)-L 

zq z for X < ^ < L 


( 87 ) 


( 88 ) 


and therefore to densities 

( 5 [l — tanh(^(^ —x + Ll^))] for^<x 
Pk = I I for k = x (89) 

5 [l — tanh—x + Ll^))] for^>x 

where £ = In^ and k = In z/E corresponding to z = q’^. These measures are closely 
related to the type-II shock measures defined in Q. On coarse-grained scale with 
^ = k/L and =xlL the macroscopic density profile p(^) has a discontinuity 
at ^ with amplitude A-* := p+ ~ p_ = tanh (£(«•+ l)/2) — tanh (£(«•+ l)/2), 

where = lim£_j.op± e). In forward (clockwise) direction and for E > 0 this 
is an upward step with fugacity ratio z^ jz^ = cf', corresponding to a shock profile 
for the ASEP (with positive bias E > 0). For strong asymmetry E = eL one has 
near k = x — L{1 + k)/2 mod L a smoothened downward “step” with an intrinsic 
width oc 1/e on lattice scale. On macroscopic scale this constitutes an antishock at 
position — (1 + k)/2 mod 1. 

For K = 2 the SAM | Px,y ) (l85l l with 1 <x<y<L has local fugacities 


l{lk-x-y-\-L) 

zq z 

tor 1 < AC < X 

00 

for k = x 

2(2k-x-y) 

zq z 

for X < k <y 

00 

for k = y 

2(2k—x—y—L) 

zq z 

tor y < k <L 


corresponding to densities 


Pk 


' \ [1 -tanh(f (2k-x-y + L^))] 

1 

< 7 [l — tanh(j(2k —X —y+ Lf))] 

1 

^ i [1 -tanh(f (2k-x-y + L^))] 


for 1 < k < X 
for k = x 
for X < k <y 
for k = y 
fovy <k <L 


(90) 


(91) 


On macroscopic scale this density profile has two shock discontinuities at = 
XIL mod 1 and mod 1 ■ Both fugacity ratios are of magnitude . For strong 

asymmetry £ = eL there are two antishocks at ^ 2 )('^ + 2)/4 mod 1 

and ^2 = + ^ 2 )/^ ~ ('^ + 2)/4 mod 1. 
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3 Results 

Before stating the new results we recall the duality relation for the ASEP with re¬ 
flecting boundary conditions derived in 13^ . Eq. (3.12). We reformulate this duality 
relation slightly and correct a sign error in Eq. (3.12) of ||33l. We also give a new 
proof, parts of which are then used to prove the new results given below. We also 
present a generalized and slightly reformulated version of the intertwiner relation 
Eq. (2.62b) of lEOl for perdiodic boundary conditions, also with a correction of 
some sign errors in that formula. 

Theorem 1. (Schiitz, 4331? ) The ASEP with reflecting boundary conditions and 
asymmetry parameter q is self-dual w.r.t. the duality function 


j=i 


where 




Remark 1. Because of particle number conservation also 


(92) 


(93) 


D(x,77) (94) 

;=i 


is a duality function with the rarticle numbers N^r]) Q and Nx{ri) (|4]i. This is the 
duality function (3.12) of IESE 

Proposition 1. Let Hf, •, •) be the conditioned generator (|29ll of the ASEP with pe¬ 
riodic boundary conditions and let rjK € Qk be any configuration with K particles. 
Then for 0 <n < L — K one has the intertwining relation 

[iS^iq,a)rH{q,a,q^"l5^)-H{q,a,l5^){S^{q,a)r]\iqK)=Q (95) 


with 

(96) 

and the generators S^{q,a) (|62]l ofUq[Q{{2)\. 

Remark!. Defining the duality matrix = lK^^n{S^{q, with the projec¬ 

tor (l45l l and using 438]) the intertwiner relation 495l l can be expressed as the duality 
relation 

D^^"HKiq,a,q^'’p^) = {HK^n{q,a-\p^)fDf^'' (97) 

with the projected generator 4461 ). We shall focus on the formulation 495] ) of this 
duality. 


^ Notice a sign error in front of the term 2kj in Eq. (3.12) of (33) and pay attention to the different 
convention q . 
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Remark 3. For a = 1 this is the result (2.62b) of |[^ FI The proof of Proposition ([Til 
is entirely analogous to the derivation given in ll^ since the generalized form (l95T l 
follows trivially from the result of ll^ for a = 1 through the similarity transfor¬ 
mation (l33l l. Some ingredients of the proof, with sign errors in ll^ corrected, are 
presented in the appendix. 

We focus now on global conditioning {a ^ q, j3 = 1) and local conditioning 
oc = q, j3 1. The main results of this work are the following theorems. 

Theorem 2. Let H§ := lNH{q,q^^^ ,1)1^ be the generator ( 1461 ) of the globally 
conditioned ASEP with N particles and periodic boundary conditions and driv¬ 
ing strength s = —2K/L\nq. Furthermore, let |/i^) = IatI/Tx) be the unnormalized 
shock-antishock measure ( liS5l ) restricted to N particles and 

(98) 

with /T = |x|. Then 

l<(0)= E ^"(y4|x,0)|O (99) 

yeiiK 

where P^{y,t\x,0) := (y is the conditioned K-particle transition proba¬ 

bility from X toy at time t with driving strength s' = —2N/L\nq. 

Remark 4. The signihcance of this result lies in the fact that the conditioned evolu¬ 
tion of an A^-particle SAM is fully determined by the conditioned transition proba¬ 
bility of only K particles, in analogy to the evolution of shocks in the inhnite lattice 
explored in 

Remarks. For /T = 1 a related result was obtained in 0 for a normalized and 
slightly different dehnition of the shock measures. The proof of 0 is by explicit 
computations relying on the presence of a single shock. The present proof for the 
generalized K>\ shows that the mathematical origin of the conditioned shock mo¬ 
tion is the duality relation ( i95] l. 

Theorem 3. Let H§ := lNH{q,q,q^^^)'^N be the generator ( 1461 ) of the locally con¬ 
ditioned ASEP with N particles and periodic boundary conditions and boundary 
driving strength s = —2K\t\q. Furthermore, let |/i^) = IatI/Ix) be the unnormal¬ 
ized shock-antishock measure (|S2l) restricted to Nparticles and 

K(0):=e-^-'|Alx^) (100) 

with fif = |x|. Then 

\fi!^it))= Y. (101) 

where .P^(y,f|x,0) := (y |e^^«^|x) is the boundary-conditioned K-particle transi¬ 
tion probability from x toy at time t with driving strength f = —2N\nq. 


5 Eqs. (2.62a) and (2.62b) of (30) have some sign errors which are corrected in Proposition o. 
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4 Proofs 


4.1 Proof of Theorem^ 


Proof. We first note 

Lemma 1. Let 


L c+ 






'Hr. 


Then for a configuration x € i2Ar w/ffi = |x| particles one has 


|x| 


(x|5= 


i=l 


( 102 ) 


(103) 


The proof is completely analogous to the proof in Oof (164) withy = 0. 

Now we observe that with the reversible measure ( fTST i and with ( l74l i we can write 


D{x,ri) = n \x){s\flQ,,\ri)=f{ri)=g'^{x). 

i=\ 

Then the following chain of equalities holds and proves the theorem: 

= i:;r-i(x)(5|na,|^)(^|e-^'|77) 

^ ,=i 

= ?r^Hx)(^inav,e^'^'|77) 

/=! 

= ;r-i(x)(x|5e-«'|t7) 

= ;r-^(x)(x|e-"'5|T7) 

= ;r-i(x)^(x|e-«'|y)(y|5|t7) 
y 

|y| 

= L (x|^^‘e-"'7f|y);r-ny)(^iriSv,|j?) 

yeQfj i=l 

yeilN 1=1 

= L ^’’(y)(y|e-"^|x) 

yeQfj 

=■■ {8Ht)r 


(104) 

(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 

( 111 ) 

( 112 ) 

(113) 

(114) 
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The following ingredients were used: Eqs. ( I105l l and (II 141 i: The expressions ( fTTT i 
for expectations; Eqs. ( 1106b and ( II 13b : The expression (1104b for the duality func¬ 
tion; Eq. (1107b : The expressions ( iTTb and the representation ( l47b of the unit matrix 
of dimension 2^; Eq. (1108b : The expression (1103b of part of the duality function 
in terms of the symmetry operator 5; Eq. (1109b : The t/^[s 1(2)]symmetry (1^ : Eq. 
(II 10b : Particle number conservation and the representation (l45T l of the unit matrix 
in the subspace of N particles; Eq. (II 1 lb : The diagonal matrix representation of the 
reversible measure (fTsT i: Eq. (II 12b : Reversibility d^Tb . □ 


4.2 Proofs of Theorems ^and^ 


Before we set out to prove Theorem|2]and Theorem|3]we show that the SAM (l85b 
can be generated by the action of the particle creation operator t/^[s[(2)]. 

Lemma 2. Let xbe a configuration ofK = |x| particles and let 

jl^ := jl^5^L^^ri{k],N’ T](yt),A (115) 

be the SAM’s defined by ( |iS2b , ( |iS5b restricted to N > K particles. Then the vector 
representations \ ft^) := IatI /ix) and \ ) := 1a:| fix) can be written as 




lAx^) = 


[N-K],: 


I^A) V- J 


N-K 


[N-K],: 


|x) 


in terms of the generators ( 162b ofU,[sl{2)] and the transformation ( 133b . 
Proof Erom Lemma [T] and (l66b one finds 


(116) 

(117) 


L-K 


E 

7?=0 


(5 {q \q))'\ , A - 

—u-i —hi-n? 


J9' 


j=i 




(118) 


Notice that for z = 1 one has j z = 1) = j s). 

The transformation ( |3^ yields V( 7 )|x) = gives 

,q) = (A). Putting this together and using (1651) turns 

(II 18b into 




/!=0 


[n]q\ 


./=i 


( 119 ) 
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Now we choose A 


2N 2{K-N) . 

q'Tr and 7 = q to obtain 


2(K-N) ^)) 

V{q^) E ^ ^ ' 

«=0 


W 9 ! 


/i: 


>=i 


( 120 ) 


Finally one applies the projector on both sides of the equation. On the l.h.s. 
this projects out the term with n = N — K, corresponding to the r.h.s. of dl 171 ). On 
the r.h.s. the projection allows us to substitute the number N in the first power of q 
by the number operator N (l34l) . Thus the terms proportional to L + 1 cancel and the 
expression (l 86 l) remains under the projection operator. Therefore the r.h.s. is equal 
to I ). Similarly one chooses 7 = A = 1 to obtain (II 16b . □ 


4.2.1 Proof of Theorem 12 

Now we are in a position to prove ( l99b . 

Consider a /T-particle configuration x and the duality relation ( l95l l with n=N — K 
and j3 = 1 : 


2^-\\\N-K 


{S {q,q^i- *)) 
[N-K],\ 


-H{q,q^T-\q^^-2^)\x) 


With the transformation (l3^ with 7 ^ = one uses drill to cast this in the form 

’ [N-K]q\ ' ' 


[N-K]q\ 


V-\y)H{q,q^T-^,l)V{y)\x) (122) 


or, alternatively, 




[N-fi:]..! 


?• 


2‘^-\\\N-K 


= V-\y) 




H{q,q^T-\\)\^). 


[N-K]q\ 

Applying (ITTI) . dCTl) . dSSl) this turns into 


(123) 
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[N-K]J- 


H{q,q^-^T^l)V{y) 


|x) 


9- 


(S-{q-\q^-^T)f-^ 




(124) 


where |x) = V^^( 7 )|x) is an arbitrary /T-particle configuration. 

Since H conserves particle number, this relation remains valid for any power of 
H. Thus we find 




[IV-/:],! 




[N-K],\ 


[N-K],\ 


(125) 

(126) 


where in the last equality we have inserted the unit operator restricted to /T-particle 
states. Using (II 171 i of Lemma|2]then proves (l99l l. □ 


4.2.2 Proof of Theorem |3] 

The proof of Theorem|3]is similar. For a = q^^ where j3 = q^^ one has 

Applying space reflection dJTl i. dfiTl i this becomes 


- H{q,q,q ^^)\x) = H{q,q,q - - - - — 

Here we dropped the tilde over the configuration x since it is arbitrary. 


|x). (128) 


Projecting on N particles and iterating this duality over powers of hJ^ yields 
[A-/:],! ' ^ i ' ^ 


[A-A],! 


Inserting the unit operator restricted to A-particle states and Using (II 16b of Lemma 
|2]then proves (llOlb . □ 


Acknowledgements GMS thanks DFG and FAPESP for financial support. 
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Appendix 

We present some details of the proof of Proposition ([T]i which are not shown in ll^ 
and from which Proposition ([T]i follows by the similarity transformation (l3^ . 

We define eL(-, •, •) := /li.i (•, •, •), see OTT i. By explizit matrix multiplications one 
finds from the relations (l25T l for the bulk operators 

S^iq,a)eL{a',q',li) = eLia\q',l5q-^)Sfiq,a) 2<k<L-l (130) 

eL(a',q\l5)S^{q,a) = S^{q,a)eL{a',q',pq^) 2<k<L-l (131) 

and for the boundary operators 

S^{q,a)eL{a’,q',l5) = [{q'y'^a+VL- a’l5vial]q^^'' (132) 

Sl{q,a)eL{a',q\l5) = [q'a{hL-{a'py^nia^]q^^' (133) 

Sl{q,a)eLia',q,p) = [q'via,^ - {a'py^a^Vilq^' (134) 

Sl{q,a)eL{,a',q' (135) 

and 

eya',q',fi)syq,a) = [q' o^il- a'fihiol]q^^'' (136) 

eya',q',^)S\{q,a) = [(gr')-i(7-f)^ _ (a'j3)-ibjo);] ^^'^(137) 

eya\q',l5)Sl{q,a) = [{q'y^fiial -{a'liy^a^hL] q^' (138) 

eya',q',p)Sl{q,a) = [^'■Oio); - a'iSafOi](139) 

Consider now q = q' and a = a'. From the quantum algebra symmetry and from 
the previous relations one obtains (omitting the q, a-dependence) 

S^H{j5)-Hil5')S^ = syyi5)-eyi5')S^ (140) 

= [eyq-^i5)-eyy)]'fsf 

k=2 

+ {Sf + 5±) eyp)-eyp') {Sf + 5±). (141) 

Observe that 

= Si =q^^^a-^-^^-^'>a{q-^\ (142) 

S+ = SZ = (143) 

and the auxiliary relations 

a+eyp) = q^^a^VL-apyay 
eyp)a^ = qayiL-a^niap 


(144) 

(145) 
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a,^eLil5) = qvial-{apy'^a^VL, (146) 

eL(j3)a+ = q-^nia+ - {apy^a^fiL, (147) 

and 

a^eyp) = qa^nL-iapy'nia[, (148) 

eL(j3)(Tf = q^^a^;vL-{al5y^Via^ (149) 

0);eL()3) = ^^'niC7^-aj3(JfnL, (150) 

eL(j3)o); = ^Oio);-ajSafOi. (151) 

Thus one obtains 


{Si+S+)eL{l5)-eyi5'){Si + S+) 

= A+{l5,y)y^^a^^^q-^"-^ +B+(j3,j3')j3-‘/^a-^/V'+‘ (152) 


with 


A+{yy) 

B+{i5,y) 


- (7+ [q 'vl- qni] - (J? 


(apy/2 

CT+ 

(a)3)'/^ ' 




1/2 


y 

qvi-q-pni 


(153) 


-1 ^ - -i.'\ 

q jjnL -q vl 


[q '« i -? 11 (} 54 ) 


and 


{s^+sz)eyi5)-eyy){si+sz) 

= A-{yy)l5-^/^a-^/^q-^"+^ +B-(j3,j3')j3‘/^a^/V'^‘ (155) 


with 


A-{p,n = 


(a)3)'/^ 


, 1/2 


(7j [qht-q ^Dl]- 


q 


1/2 


q ^ni-q '^'0i(i56) 






1/2 -1 ^qjVL-qni 


(aj3)'/^ ^ 

fll/2 

(157) 


With the choice P' = q (1141b reduces to 
S^HiP) -H{q-^p )5± = (5f + Sp) eyP)- cl[ q-^P ) + 5±). 


( 158 ) 


For 5+ the r.h.s. reduces to 
7 I /2 


q 


(apy/^ 


(^ 1 + [q 'vl- qni] - “ q ‘«l] | 




L/2 -S--1 


_^-l/2^-L/2^y=+l 


X 
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With (I142l i. ( 11431 ) one thus arrives at 

S+H{j5)-H{q-^P)S+ = [q-^VL-qnL\Si [l - 
+ [qV\~q^^hi\S^ 1 - 


(159) 


Notice that the action of the pseudo commutator on states with particle number 
satisfying 




,L-2N+2 


= l5a^ 


(160) 


vanishes. 

Similarly one obtains for the r.h.s. of 


- r - 1 ' - 1 9''^^ - r ^ -i^ i 

[q VL-qnL]-j^^^cy^[qVi-q uj] 




1/2 


l^l/2^L/2^~-l _ i^-l/2^-L/2^-S’~+l 


a ^'^q 


which yields 


S-H{j5)-H{q-^j5)S- = -[q-^VL-qfiL]SY l-pa^q 


L2S~-2 


— [qVi—q *ni] 5 ^ ^Oi ^q 


l^-L„-2S=+2 


■ (161) 


Notice that the action of the pseudo commutator on states with particle number 
satisfying 


q 


-L+2N+2 


= l5a^ 


(162) 


vanishes. 

In compact form (1158b can thus be written 

S^H{p)-H{q^^j5)S'^ = ±[q^^VL-qnL\Sf 1 - 
±\q\)l—q^^hi\Sf^ 

One can iterate. E.g. for (5^)^ one obtains 

= {\+q-^)[qhL-q-'^VL]S^ [l 


. (163) 


L-l 

+ (\+q^^)[q^'^ni-qVi\ ( ^ 

k=2 


\-^n,-L„- 2 S'+A 


■ (164) 


Iterating further as in ll^ one arrives at Proposition[T] 
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